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AERONAUTIC SYMBOLS 


1. FUNDAMENTAL AND DERIVED UNITS 



Symbol 

Metric 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbrevia- 

tion 


l 

meter 

m 

foot (or mile) 

ft (or mi) 


t 

second— 

s 

second (or hour) 


Force- - - . 

F 

weight of 1 kilogram 

kg 

weight of 1 pound- 

lb 

Power 

P 

horsepower (metric) 


horsepower 

hp 


V 

/kilometers per hour— . . 

kph 

milesper hour 

mph 



(meters per second 

mps 

feet per second. 

fps 


”| • 2. GENERAL SYMBOLS 

W Weight=mgr p - Kinematic viscosity 

g Standard acceleration of gravity=9.80665 m/s 3 p Density (mass per unit volume) 

i or 32.1740 ft/sec 2 Standard density of dry air, 0.12497 kg-m -t -s 3 at 15° G 

v , W and 760 mm; or 0.002378 lb-ft -4 sec 2 

m Mass— ^ Specific weight of “standard” air, 1.2255 kg/m® or 

I Moment of inertia=mi*. (Indicate axis of 0.07651 lb/eu ft 
; radius of gyration k by proper subscript.) 
a Coefficient of viscosity 

j 3. AERODYNAMIC SYMBOLS 


s 

s„ 

G 

b 

G. 

A 

V 

4 

L 

D 

D * 

Dt 

D , 

G 


Area 

Area of wing 
Gap 
Span 
Chord 


b a 


Aspect ratio, g- 
True air speed 
Dynamic pressure, ^pF* 

Lift, absolute coefficient 


Drag, absolute coefficient C D — 


D 


A 


Profile drag, absolute coefficient 


D t 


Induced drag, absolute coefficient A>,=^ 
Parasite drag, absolute coefficient C Dp =^ 


C 


it 

Q 

Q 

B 


a 

e 

a„ 

a« 


Gross-wind force, absolute coefficient Cc=^g 


Angle of setting of wings (relative to thrust line) 
Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 


Reynolds number, p-^ where l is a linear dimen- 


sion (e.g., for an airfoil of 1.0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 

Flight-path angle 
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THEORETICAL MOTIONS OF HYDROFOIL SYSTEMS 


By Frederick H. Imlay 


SUMMARY 

Results are presented of an investigation that has been 
undertaken t'o develop theoretical methods of treating the mo- 
tions of hydrofoil systems and to determine some of the important 
parameters . Variations of parameters include three distribu- 
tions of area between the hydrofoils, two rates of change of 
down, wash angle with angle of attack, three depths of immersion, 
two dihedral angles, two rates of change of lift with immersion, 
three longitudinal hydrofoil spacings, two radii of gyration in 
pitching, and various horizontal and vertical locations of the 
center of gravity. Graphs are presented to show locations of the 
center of gravity for stable motion, values of the stability roots, 
and motions following the sudden application of a vertical force 
or a pitching moment to the hydrofoil system for numerous sets 
of values of the parameter's. 

The lateral stability of tandem-hydrofoil systems is briefly 
discussed, and values of the lateral stability roots are presented 
for two longitudinal hydrofoil spacings and two vertical loca- 
tions of the center of gravity. 

The analysis indicates that if only the longitudinal motions of 
a hydrofoil system are of interest the present theory should pro- 
vide satisfactory predictions. An adequate theory for the 
lateral motions, however, must treat the longitudinal and 
lateral motions in combination. The conclusions based on the 
investigation are that a large longitudinal spacing between the 
hydrofoils, a large rate of change of lift with depth of immersion, 
and, a horizontal location of the center of gravity near the center 
of the region of stable locations are important contributions in 
the, attainment of desirable characteristics for the longitudinal 
motion. An appendix gives an outline of the methods of 
theoretical treatment used and presents methods used- in com- 
puting the required stability derivatives. 

INTRODUCTION 

The use of hydrofoils as an alternative to planing bottoms 
or hulls for the support of craft operating on the surface 
of water has been of interest for some time. (See reference 1.) 
Guidoni advocated the use of hydrofoils as a means of 
improving the take-off and rough-water performance of sea- 
planes as early as 1911. (See reference 2.) Some of the advan- 
tages claimed for hydrofoils over planing bottoms are a better 
ratio of lift to drag on the water and less sensitivity- to 
irregularities of the water surface. In addition, if hydrofoils 
are used, the hull lines can be designed to favor good aero- 
dynamic rather than good hydrodynamic characteristics, and 
by retracting the hydrofoils the aerodynamic performance 
can be even further improved. In spite of the evident 
advantages of these devices and the attention that they 
have received, no published work is known to exist on the 
stability of motion for systems employing hydrofoils. 


The present paper deals theoretically with the behavior 
of a system supported solely by hydrofoils and is a first 
approach to the problem of developing methods of theoretical 
treatment for the more general case where the interaction 
of hydrofoils, hull, and aerodynamic surfaces have to be 
taken into account. The treatment is based on the theory 
of small oscillations and involves assumptions customarily 
made in applying the theory. (See reference 3.) 

Definitions of all symbols used are listed at the beginning 
of the appendix. 

LONGITUDINAL MOTIONS 

The longitudinal motions of a number of hypothetical 
hydrofoil systems were investigated by means of calculations 
based on the theoretical treatment presented in the appendix. 
All the computations were for systems composed of two 
similar hydrofoils of rectangular plan form and rectangular 
tips. The hydrofoils were arranged in tandem and had an 
aspect ratio of 6 and a total hydrofoil area of 0.188 square 
foot. (See. fig. 1.) The systems were assumed to have 
a mass of 0.25G slug and to operate at a velocity of 20 feet 
per second in water having a density of 1.97 slugs per 
cubic foot. The mass of the system was assumed to include 
all items such as structure and additional mass effect. For 
systems with dihedral the hydrofoil area, aspect ratio, and 
span were based on the part of the hydrofoil immersed 
during the initial undisturbed motion, although unwetted 
parts of the hydrofoils were assumed to project above the 
water far enough to ensure that the tips were never immersed 
during disturbed motions. (See fig. 2.) Most of these 
dimensional characteristics of the hydrofoil systems were 
chosen to facilitate comparison of the theoretical motions 
with the results of contemplated experimental tests. Changes 
in the other parameters were made to determine their effects 
on the stable regions, the stability roots, and the motions 
resulting from disturbances. 

EFFECT OF PARAMETERS ON STABLE REGIONS 

The stable region, as used in the present paper, indicates 
permissible locations of the center of gravity relative to the 
hydrofoils if the longitudinal motions are to be stable. The 
stable region alone, however, gives no quantitative indication 
of the degree of stability. The stable region is bounded by 
lines that are the loci of center-of-gravity locations for which 
neutrally stable longitudinal motions occur. The positions 
of the boundary lines, and hence the size of the region, vary 
with changes in the parameters of the hydrofoil system and 
thus suggest variations of the parameters that may be of 
practical interest for more detailed study. 
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(a) Arrangement 1. 

(b) Arrangement 2. 

(e) Arrangement 3. 

Figure 1. — Hydrofoil plan -form arrangements assumed. 

The type of unstable motion occurring just outside the 
boundaries has been noted for each of the stable regions in 
figures 3 to 9; thus, for each stable region, center-of-gravity 
locations beyond the rear boundary lead to an unstable 
divergence, and in most cases unstable oscillations occur for 
locations beyond the front boundary. The rear boundary 
is always located farther to the rear of the front hydrofoil 
than would be the case for a similar pair of airfoils because 
of the additional damping introduced as a result of the 
sensitivity of the hydrofoils to depth of immersion. 

In addition to the selection of a center-of-gravity location 
that lies within the stable region in order to meet the require- 
ments for stability, certain supplementary practical factors 
must be considered. For example, negative lift on either 
hydrofoil should be avoided; otherwise momentary uncover- 
ing of the hydrofoil (as by a wave trough) will be followed 
by nosing-over if the rear hydrofoil is operating at negative 
lift, or nosing-up if the front hydrofoil is operating at 




Figure 2. — Definition of symbols for a representative hydrofoil system. 

negative lift. Furthermore, the longitudinal location of the 
center of gravity is also restricted by the maximum positive 
lift obtainable, and may be influenced by the desirability of 
operating the hydrofoils near their maximum lift-to-drag 
ratios. The net effect of such restrictions is to reduce the 
usable part of some of the computed stable regions shown 
in figures 3 to 9. 

In the present study, where the effects of power are 
neglected, the vertical center-of-gravity location selected ap- 
pears to be of secondary importance, low locations being 
somewhat advantageous. The effects of power, however, 
will undoubtedly have an important bearing on the choice 
of the vertical centcr-of-gravity location. 

Distribution of area. — The effect of the distribution of area 
between the two hydrofoils on the extent of the stable region 
is shown in figure 3. The plan-form arrangements assumed 
for the three distributions treated are shown in figure 1. In 
arrangement 1 the hydrofoils were identical; in the other 
two arrangements the ratio of the distribution of area was 
1 : 4 and the arrangements differed only in the location of the 
larger hydrofoil. All the arrangements had the same total 
hydrofoil area of 0.188 square foot. The horizontal distance 
between the assumed hydrodynamic centers of the hydrofoils 
for all arrangements was lO.Oci, where c i is the chord for the 
arrangement with two equal hydrofoils, and the assumed 
hydrodynamic center was located at the quarter-chord point 
of the center section. All the hydrofoils were assumed to 
be immersed l.Oci at the hydrodynamic center during the 
initial undisturbed motion. 

Figure 3 shows that the configura tion with the small surface 
ahead (arrangement 3) gave the largest useful stable region. 
The rearward extent of the stable region for the arrangement 
with two hydrofoils of equal area (arrangement 1) was con- 
sidered adequate, however, and because this arrangement 
permitted certain simplifications in the calculations, it was 
used for the rest of the work. The configuration having the 
main surface ahead (arrangement 2) would, from theoretical 
considerations, be the most efficient arrangement for de- 
veloping lift but has a considerably more limited range of 
stable center-of-gravity location than do the other arrange- 
ments. 
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Figure 3— Stable regions for three area distributions. r = 0°; e=0; / = l().0c>; z 0 =l.OOa; ATr=6.67ci. (For this figure c\ is the chord for arrangement l,in which the two hydrofoils are identical.) 



Figure 4.— Stable regions for two downwash angles. r=0°; Z=l0.0ci; z„=1.00ft; IiTr= 6.67c i, 
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Figure 5.— Stable regions for throe immersion depths. r=0°; Si-Si; e=0; Z=10.0ci; 2£V=6.67ci. 



Figure 6.— Stable regions for two dihedral angles. Ss=S\; €=0; Z=10.0ct; z o =1.00ci for r=0°, and 1.74ci for r=30°; -KV=6.67ei. 
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Figure 8. — Stable regions for two longitudinal spacings. r=30°; Sa^Si] e=0; z 0 =1.74ci; Ky=S.Q7ci. 
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Figure 9. — Range of stable ccntcr-of-gravity location for two values of Ky. r=30°; Si~S\, e=0; i — lO.Oui ; z<,=1.74ci. 


Rate of change of downwash. — In a tandem-hydrofoil 
system, the downward velocity produced in the fluid by the 
front hydrofoil reduces the effective angle of attack of the 
rear hydrofoil by the amount, of the downwash angle e. The 
downwash angle is a function of the lift on the front hydrofoil 
and hence varies with angle of attack. The rate of change 
of downwash angle with angle of attack, which is the factor 
of interest from the standpoint of stability, will be repre- 
sented by the symbol The value of e a will probably be 
intermediate between zero and the theoretical ultimate maxi- 

2 d(C ) 

mum ^ * but to determine the value accurately 

would require an investigation of downwash near a free 
surface. Corresponding limiting values of e, which are given 
instead of « a in the figures for the sake of brevity, are zero 
and twice the induced angle of attack a,-. In order to show 
the influence of the rate of change of downwash on the nature 
of the stable region, computations were made for these two 
extremes, and the results for a system having two equal hy- 
drofoils are shown in figure 4. An increase in the variation 
of downwash with a. shifts both boundaries forward without 
appreciably altering the size of the stable region. 

The effect of downwash for the other hydrofoil arrange- 
ments was found to be similar to that indicated by figure 4 
for the arrangement with two equal hydrofoils. Because 
there was no pronounced change in the size of the stable re- 
gion with change in downwash, the condition of zero rate of 
change of downwash with a was assumed in most of the 
remaining calculations. 

The true boundaries of the stable region for the system 


treated in figure 4 lie some.vhere within the bands defined 
by the boundaries for e=0 and e=2a,-,but accurate definition 
of the boundaries requires that e be known. Conservative 
estimates will be obtained, when the value of e is not known, 
if the assumptions are made that e—2a t for computing the 
location of the rear boundary and that e=0 for the front 
boundary. 

Depth of immersion. — The lift and drag obtained from a 
hydrofoil depend upon the depth of immersion z 0 of the 
hydrofoil in the water. Because appreciable change in the 
depth, of immersion may occur under normal operating condi- 
tions, computations of the stable regionjwere made for 
immersion depths of 0.5ci, l.Oci, and 1.5ci. (See fig. 5.) 
Limits of the stable region were not altered to any important 
extent by the assumed changes in the depth at which the 
hydrofoils operate. 

Dihedral angle. — The effect on the stable region of increas- 
ing the dihedral angle r of the hydrofoils from 0° to 30° is 
shown in figure 6. Both the front and the rear boundaries 
of the stable region were affected by the dihedral in such a 
way that the increase in dihedral increased the size of the 
stable region. 

Increasing the dihedral angle from 0° to 30° resulted in 
an associated increase in vertical damping. It appeared 
reasonable that the improved stability obtained by changing 
the dihedral might have resulted from this increased vertical 
damping; consequently the effect of arbitrarily increasing 
the vertical damping for the hydrofoils with a dihedral angle 
of 0° was studied and the results are discussed in the next 
section. 
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Rate of change of lift with immersion. — If the depth of 
immersion of a hydrofoil is changing, the lift is also changing, 
and the rate of change can be expressed by the vertical- 
damping derivative bC L /bz'. It is believed that the increased 
stability which accompanied the increase in dihedral angle 
from 0° to 30° (discussed in the preceding section) may 
have been brought about by the resulting increase in the 
value of dC L /dz'. Inasmuch as a further increase in dihe- 
dral angle would decrease the value of the derivative, an 
explanation of the increase in dCJdz' when the dihedral 
was changed from 0° to 30° may be of interest. 

In order to avoid the mathematical difficulties of treating 
discontinuous derivatives the assumption was made in the 
present study, for the case of hydrofoils with dihedral, 
that a normally inactive part of the hydrofoil extended 
sufficiently far above the water surface to keep the hydrofoil 
from being completely immersed at any time during dis- 
turbed motion. (See fig. 2.) As a result of this assumption, 
hydrofoils with dihedral have a larger variation of lift with 
change in depth of immersion than do hydrofoils with 0° 
dihedral because of the increased area brought into action 
when the hydrofoil sinks deeper into the water. This 
variation in active area becomes greater as the dihedral 
angle becomes smaller. 

The effect on the stable region of increasing the value of 
dCJdz' for each hydrofoil to twice the value that the hydro- 
foils had with 30° dihedral, but having other characteristics 
the same as for 0° dihedral, may be seen by comparing figures 
6 and 7. Doubling the value of dCV./ds' shifts the rear 
boundary of the stable region back considerably and produces 
pronounced changes in the front boundary. The former 
boundary for unstable oscillations now becomes an unstable 
“hump” in the region with a new front boundary ahead of 
the hump. The new forward boundary represents conditions 
for an unstable divergence, but the boundary is too far ahead 
of the front hydrofoil to be of any practical interest. 

Longitudinal hydrofoil spacing. — The effect on the stable 
region of increasing the longitudinal spacing of the hydro- 
foils from 10ci to 20ci is shown in figure 8. The larger spacing 
results in a very large increase in the stable region and in the 
replacement of the front boundary that indicated unstable 
oscillations by a new front divergent boundary. The new 
front boundary is well ahead of the front hydrofoil, which is 
the practical limit of forward center-of-gravity location. 

The absence of a boundary for oscillatory instability for 
the system with a spacing of 20c! suggests that the large 
amount of damping in pitching for this spacing, relative to 
the pitching radius of gyration K Y , might result in over- 
damping and thus prevent the system from having any oscil- 
latory motion. Calculations with K Y reduced to give a 
similar relation between inertia and damping for the small 
spacing of 10ci, made to check the hypothesis, showed that 
oscillations were still obtained; thus, it appears that the 
absence of unstable oscillations for the larger spacing does not 
signify inability of the system to have transient oscillations. 

The pronounced increase in the size of the stable region 
when the longitudinal spacing of the hydrofoils is increased 
indicates that a large spacing is desirable in order to minimize 
the effects of unavoidable changes in center-of-gravity 
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location encountered in practice. In a previous section entitled 
“Distribution of area,” a spacing of 10ci was used in the 
calculations made to study the effects of distribution of area. 
If a larger spacing had been used, it would possibly have 
resulted in a sufficient gain in the size of the stable region for 
the arrangement with the large hydrofoil forward to make 
this configuration of practical value. 

Radius of gyration in pitching. — The marked increase in 
permissible horizontal center-of-gravity movement when 
K y is reduced is indicated in figure 9, where the stable range 
of horizontal center-of-gravity location is shown for zero 
vertical elevation of the center of gravity with K y reduced 
to one-fourth the value used previously. The pronounced 
effects of reducing K Y indicate that increased values of K y , 
which are more likely to be used, should receive attention 
because of possible adverse effects on the characteristics 
of the longitudinal motions. 

EFFECT OF PARAMETERS ON STABILITY ROOTS 

When the equations of motion are solved, the motion is 
obtained as the sum of a series of components called modes. 
Stability roots, which indicate the degree of stability of the 
various modes, can also be obtained from the equations of 
motion without effecting a complete solution of the equa- 
tions. A more detailed discussion of the significance of the 
stability roots is contained in the appendix of reference 4. 
Information obtained from the stability roots is most useful 
when the relative magnitude or importance of the various 
modes is known, because the roots then provide a clue to the 
nature of the complete motion. 

In the present analysis, four stability roots X arc obtained 
from the longitudinal equations of motion and are distin- 
guished by the subscripts 1 to 4. The nature of the roots 
changes with variations in the parameters of the hydrofoil 
system. A typical variation in the real parts of the roots 
is shown in figure 10 (a). In general, when the magnitudes 
of any two real roots become equal, the two real roots arc 
replaced by a conjugate pair of complex roots, each having 
the same magnitude for the real part. Thus, such pairs of 
complex roots in figure 10 (a) are indicated by a double line 
and an appropi’iate modification of the subscript. The 
magnitude of the real part for such complex pairs of roots 
should be read off the plots at the center of the double line. 

For every real root obtained from the equations of motion 
the complete solution will contain an aperiodic mode, or 
component, of the motion. Likewise, for every pair of 
complex roots the motion will contain an oscillatory com- 
ponent. When the magnitude of the real part of any of 
the roots passes through zero, the motion becomes unstable. 

Horizontal center-of-gravity location. — The effect of chang- 
ing the horizontal location of the center of gravity on the 
real parts of the stability roots is shown in figure 10 (a) for 
a system of two equal hydrofoils with 30° dihedral. For 
center-of-gravity locations ahead of the hydrodynamic cen- 
ter of the front hydrofoil, two real roots Xi and X 2 and a pair 
of complex roots X 3 . 4 exist. When the center of gravity is 
2.16cj ahead of the front hydrofoil the X 3 ,i roots are unstable, 
which indicates that the center of gravity has reached the 
forward boundary of the stable region. As the center of 
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gravity is moved rearward, the stability slowly improves for 
the oscillatory component of the motion represented by the 
X3.4 roots. Meanwhile the magnitudes of the Ai and X 2 
roots approach each other and become equal when the center 
of gravity is about 1.5ci behind the front hydrofoil. With 
farther rearward movement of the center of gravity the roots 
are coupled as two oscillations represented by X 4 . 2 and X 3 . 4 . 
When the center of gravity is moved back to the vicinity of 
4.5ci behind the front hydrofoil rather rapid changes in 
coupling occur, which finally result in a real root Xi with a 
large amount of damping, a complex pair X 2 . 3 with moderate 
damping, and a real root X 4 with slight damping. When 
the center of gravity is moved back to a point 5.49ci belli. id 
the front hydrofoil, the magnitude of the X 4 root becomes 
zero and the rear boundary of the stable region has been 
reached . 

The behavior of the roots as the horizontal location of the 
center of gravity is changed indicates that the type of mo- 
tion caused by disturbances will be considerably influenced 
by the longitudinal location of the center of gravity. 

Rate of change of downwash. — The effect on the stability 
roots of assuming the downwash angle « to be 2a* instead 
of zero can be seen from a comparison of figures 10 (a) and 
10 (b). No pronounced change in the roots occurred with 
variation in e a , except for a shift of the pattern of root 
couplings with respect to the horizontal centor-of-gravity 
location; this result is consistent with indications obtained 
earlier from a study of the influence of e a on the stable 


region. Hence, for the rest of the work the value of e„ was 
assumed to be zero. 

Dihedral. — The influence on the stability roots of changing 
the dihedral angle from 30° to 0° is evident when figure 
10 (a) is compared with figure 10 (c). The difference in the 
rate at which the X 2 . 3 oscillation develops with rearward 
center-of-gravity movement for the two dihedral angles 
accounts for the different appearance of the right side of the 
diagram in the two figures. The most important feature dis- 
closed by the comparison is the improvement, brought about 
by the use of dihedral, in damping of the component of mo- 
tion involving the root X 4 or the complex pair X 3 . 4 . 

Vertical center-of-gravity location. — Figures 10 (d) and 
10 (e) together with figure 10 (a) show the effect on the 
stability roots of varying the vertical center-of-gravity loca- 
tion from a point on a level with the hydrofoils to a point 
lOcj above the hydrofoils. As had been indicated by the 
diagrams of the stable regions, no pronounced changes occur 
in the nature of the roots when the vertical center-of-gravity 
location is shifted. 

Rate of change of lift with immersion. — The effect on the 
stability roots of making the value of bCJbz' twice that for 
30° dihedral is evident if figure 10 (f) is contrasted to figure 
10 (a). Doubling the vertical-damping derivative caused 
marked improvement in the X 3 . 4 oscillation, which suggests 
that the similar improvement in damping obtained by in- 
creasing the dihedral angle from 0° to 30° was a result of the 
associated increase in the value of bCJbz'. 



(a) r=30°; S 2 =Si-, t=0; i=10.0c,; z„=1.74ci; JETr— 0.07c,; 2 , = 5.00n. (b) r=30°; Ss=Si; e=2 a t ; Z=10.0ci; 2 „= 1.74ci; KV=G.07d; zi=5.00ci. 

Figure 10. — Variation of real parts of stability roots with x\. 
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EFFECT OF PARAMETERS ON INDICIAL RESPONSES 

An indicia] response is the motion resulting from a unit 
force or moment suddenly applied to the hydrofoil system 
at zero time and held constant thereafter. The indicial 
responses are of interest because they are of the same general 
character as the motions produced by types of disturbance 
that are likely to be encountered in practice. 

The longitudinal equations of motion (equations (9)) in- 
volve three variables; hence three indicial responses are 
necessary to define the motion caused by any specific unit 
disturbance. The three indicial responses may be con- 
veniently represented by the symbols a z , z' z , and 6 Z for 
the change in angle of attack, vertical position, and angle of 
pitch, respectively, when the motion is caused by the sudden 
application of a unit C z force to the hydrofoil system. 
Similarly a m , z' m , and 6 m are the response factors for a 
sudden unit C m disturbance. 

The indicial responses are functions of nondimensional 
time s c , typical variations of which are shown in figure 11. 
The magnitude of disturbances actually encountered, when 
expressed in coefficient form, will usually be considerably 
less than unity; consequently, the actual motions experienced 
null be of proportionately srhaller magnitude than the indicial 



Figure 11— Indicial responses for unit Cz disturbance. r=0°; Si—S\\ <=0; Z=10.0a; 
zo=].00ci; iiTr= 6.67a; zj=5.00a; zt=-i.25a, 3.50a, or 4.80a. 


responses presented but null have the same type of variation 
with time. Values of the indicial responses after the 
disturbance has been absorbed by the system and new steady- 
state equilibrium conditions have been reached are repre- 
sented by short horizontal lines at the right side of the plots. 
Such steady-state values are not only new equilibrium con- 
ditions for sudden disturbances but also represent new trim 
conditions after gradual changes in the load condition, such 
as would result from the use of fuel. 

Horizontal center-of-gravity location. — Indicial responses 
for a unit Cz disturbance applied to a system of two equal 
hydrofoils with zero dihedral are plotted against nondimen- 
sional time in figure 11 for several horizontal locations of 
the center of gravity. Values of 'j\ used in figure 11 were 
selected to give center-of-gravity locations covering all the 
types of root coupling shown in figure 10 (c). If the center- 
of-gravity locations used in figure 11 are noted on the 
diagram of the corresponding stable region (see fig. 3) , the 
following points are evident: 

(1) A center-of-gravity location near the front boundary 
of the stable region is conducive to motions characterized 
by pronounced oscillations. 

(2) A more rearward location of the center of gravity 
reduces the prominence of the oscillations but increases the 
ultimate deviation from the attitude that existed before the 
disturbance. 

(3) For center-of-gravity locations near the rear bound- 
ary, no discernible oscillation is noted, but very large depar- 
tures from the initial condition occur. 

Comparison of the maximum deviations for the three 
center-of-gravity locations of figure 11 shows that, during 
the interval of time covered by the curves, the smallest 
amplitude of motion of the hydrofoil system occurs for the 
case with the center of gravity back 35 percent of the dis- 
tance l between the two hydrofoils. The deviation caused 
by a given disturbance rapidly becomes greater as the center 
of gravity is moved back of the optimum location, with the 
result that for such rearward locations a very slight disturb- 
ance would bring the hydrofoils to the surface or cause them 
to sink very deep into the water. Location of the center of 
gravity any appreciable distance ahead of the optimum 
location appears undesirable because of the pronounced 
oscillatory motions involved. Such motions would be both 
uncomfortable and difficult to control. 

Indicial responses for a unit G m disturbance, for the same 
conditions as for figure 11, are plotted in figure 12. The 
discussion of the effect of change in horizontal center-of- 
gravity location on the indicial responses for a unit C z dis- 
turbance also applies for a unit C m disturbance, with the 
exception that the amplitudes of the motions are least for 
the most forward center-of-gravity location considered, in- 
stead of for the middle location. The oscillations are much 
more persistent, however, for the forward location than for 
the middle location. 
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Figure 12. — Indicia] responses for unit C m disturbance. r=0°; Si=Si; «=0; Z=10.0cj; 
2®=1.00ci; JKr=6.67cj; zt=5.00ci; /i = — 1.25ci, 3.50ci, or 4.80ci. 


Figure 13.' — Indicial responses for unit Cz disturbance. r=30°; <S 2 =<Si; €—0; Z=I0.0cj; 
z 0 =1.74ci; jfiTi'=6.67ci; Zi = 5.00ci; Ji =— 2.00o, 3.50ci, 4.50ci, or 5.25ci. 


Because of the large response factors involved for either 
type of disturbance, even when the best center-of-gravity 
location is selected, motions for hydrofoils with no dihedral 
will involve large amplitudes whenever a slight disturbance 
is encountered; hence, it appears evident that such a type 
of hydrofoil will not give satisfactory performance. This 
conclusion applies only to the arrangement investigated, 
where the hydrofoils always remain completely submerged; 
and it should not be extended to cover ladder arrangements, 
for which a change in effective area with immersion depth 
produces effects similar to those for partly immersed hydro- 
foils with dihedral. 

Diehedral angle. — The effect on the indicial responses of 
increasing the dihedral angle from 0° to 30° may be obtained 
by a comparison of figures 13 and 11 for a unit C z disturb- 
ance, and of figures 14 and 12 for a unit C m disturbance. 
The figures indicate that the effect on the nature of the 
motions of changing the horizontal center-of-gravity location 
is much the same as that indicated in the preceding parts of 
the present paper. Thus, the most desirable center-of- 
gravity location appears to be about 3.50ci back of the front 
hydrofoil, as in the case for 0° dihedral angle. At any 
844045—59 3 


particular horizontal location of the center of gravity, the 
increase in dihedral causes an appreciable reduction in the 
indicial responses. The reduced sensitivity to disturbances 
when the dihedral angle was increased from 0° to 30° may 
have been a result of the corresponding increase in vertical 
damping. In such a case, as mentioned in the discussion 
of stable regions, a further increase in dihedral would have 
an effect opposite to that caused by this initial increase in 
dihedral. 

Rate of change of lift with immersion. — The effect of 
varying the rate of change of lift with immersion on the 
indicial responses for a unit C z disturbance may be seen from 
a comparison of figures 11, 13, and 15. Figures 11 and 13 
give the indicial responses for hydrofoils with dihedral angles 
of 0° and 30°, respectively; whereas for figure 15 the rate of 
change of lift with immersion is assumed to have a value 
twice that for hydrofoils of 30° dihedral angle but to have 
other hydrofoil characteristics the same as for 0° dihedral 
angle. If the case for the center of gravity at 3.50ci is 
selected in each of the figures, comparison shows the direct 
relation between good riding characteristics and a large 
value of dC' L /i)z'. It appears, therefore, that a large value 
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Figure 14. — Indicial responses for unit Cm disturbance. r=30°; Si=Si: e=0; J=10.0ci; 

Z 0 =l.74ci] JTr=6.67Cj; 0i = 5.OOa; Xi=— 2.00ri, 3.50a, 4.50a, or 5.25a- 

of dC L /c>z' should bo attained by the use of arrangements 
such as hydrofoils with dihedral for which the effective area 
changes with immersion depth, or by the use of some device 
that changes the angle of attack when the height varies. 

'Figure 16 gives data corresponding to the data of figure 15 
but with a unit C m disturbance assumed. Results for the 
several center-of-gravity locations assumed in figures 15 and 
16 indicate the same influence of horizontal center-of-gravity 
location on the motions as has been shown by the computa- 
tions summarized in figure 12. 

longitudinal hydrofoil spacing. — Indicial responses for 
either a unit C z disturbance or a unit C m disturbance applied 
to a system of two equal hydrofoils spaced 20cj are given in 
figure 17 (a). The horizontal center-of-gravity location in 
figure 17 (a) is at 0.35 /, which is the same percentage of I 
that was used in figures 13 and 14, and other conditions are 
also the same as for figures 13 and 14. Figure 17 (b) gives 
data similar to the data of figure 17 (a) except that the 
spacing has been increased to 100c!. Comparison of figures 
13, 14, and 17 indicates that increasing the hydrofoil spacing 
tends to increase the restraint in pitching and thus reduces 
the response in all degrees of freedom for pitching-moment 
disturbances, and in all but vertical motions for Z-force 
disturbances. The effect of increasing the hydrofoil spacing 



Figure 15.— Indicial responses for unit Cz /disturbance. dCiJdz' double that, for r=30° 
Si— Si] e=0; 10.0a; ifr=0.07a; 2i=f>.00a; xi = -2.00a, 0.80a, 3.50a, or G.25a. 


on the motions suggests that the spacing should be as large 
as is practical in order to reduce the response to a given dis- 
turbance. Figure IS shows the significance of 10ci, 20c!, 
and 100ci spacings if the hydrofoil systems were attached 
to a typical flying boat. 

LATERAL MOTIONS 

Lateral stability for flying boats has not generally been a 
serious problem up to the present time; hence the present 
investigation of the lateral characteristics of hydrofoils was 
brief and made chiefly to check the lateral stability of typical 
hydrofoil arrangements assumed in much of the study of 
longitudinal stability. 

In the present investigation all the lateral-stability cal- 
culations were made for a hydrofoil system consisting of two 
identical hydrofoils of rectangular plan form, each having 
rectangular tips, 30° dihedral, and an aspect ratio of 6. 
The center of gravity was assumed to have a horizontal 
location 0.35/ behind the hydrodynamic center of the front 
hydrofoil. The rate of change of downwash at the rear 
hydrofoil was assumed to be zero. The mass of the hydro- 
foil system was the same as that assumed for the investiga- 
tion of longitudinal stability. The study was confined to 
what was considered the idealized case, where the supporting 
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Figure Ifi. — Indicia] responses for unit Cm disturbance. OClIOz ' double that for l'=30°; 
52 = 5i ; e=0; z=10.0ci; .KV=6.67ci; zi=5.00ci; Xi = -2.00ci, 0.80ci, 3.50ci, or 6.25c t . 


struts have no influence on the characteristics of the hydro- 
foil system. The method of treatment for the lateral 
motions was similar to that used for the longitudinal motions 
and is described in detail in the appendix. 

The effects of changes in the vertical location of the center 
of gravity and changes in the longitudinal spacing of the 
hydrofoils on the lateral stability roots are indicated by the 
data of the following table: 


Z\ 

(chords) 

l 

(chords) 

Lateral stability roots 

2.5 

10 

0 

-0. 544=tl. 122* 

-0. 472zfc0. 204/ 

5.0 

10 

0 

-1.715=t0. 523) 

-0. 242±0. 101/ 

5.0 

20 

0 

-2. 274±1. 958* 

-0. 221 

-0. 292 


The zero root that is listed for each set of values of Z\ 
and l in the table results because the system is insensitive 
to heading; that is, the performance does not depend on the 
initial direction of travel. The remaining roots listed are 
either negative or have negative real parts in the case of 
complex roots, which indicates that all the systems investi- 
gated were laterally stable. Instability was expected in the 
two cases with the higher center-of-eravitv location, but 


apparently the stabilizing effect of the rolling moment that 
is developed when the system is banked (defined by the 
value of the derivative 2>Cyc)<f>) outweighs the effect of the 
higher center-of-gravity location. Check calculations made 
with bOi/b<j> reduced to nearly zero but with other conditions 
the same as for the second case in the table showed pro- 
nounced lateral instability. From the foregoing results the 
value of 'dC ( /c)<£ appears to have an important influence on 
lateral stability. The value of this derivative is likely to 
depend on the depth of immersion of the hydrofoils; there- 
fore it may impose a coupling between the longitudinal and 
the lateral motions and thus prevent reliable predictions of 
the lateral behavior when the longitudinal motion is ignored. 
In contrast, none of the longitudinal derivatives appears to 
be appreciably affected by lateral motions. 

The data given in the table indicate that raising the center 
of gravity and increasing the longitudinal spacing of the 
hydrofoils increase the total damping in the hydrofoil 
system, but the practical value of the increase in damping 
cannot be determined except from a study of the response 
factors involved. Such a study does not seem feasible until 
experimental checks are made on the validity of certain of 
the assumptions made in developing the theory for lateral 
motions. 


SUGGESTIONS FOR FUTURE RESEARCH 

The present study is based on the assumption of small 
displacements. Because of the nonlinearity of many of the 
derivatives involved, any appreciable departures from the 
assumed speed, depth of immersion, and other factors may 
cause marked changes in the dynamic characteristics of the 
system. Studies of maneuvers, such as take-offs, of hydro- 
foil systems may consequently require step-by-step treat- 
ment. The development of methods of studying the 
combined motions and determination of the effects of changes 
in forward speed, hydrofoil loading, and moments of inertia 
on the motions also appears desirable. For seaplanes the 
interaction of hydrofoils, hull, and aerodynamic surfaces 
must be considered. Other factors that should receive 
attention are the influence of the hydrofoil supports (partic- 
ularly on lateral motion), the effects of power, and the nature 
of the downwash near a free surface. 

CONCLUDING REMARKS 

A theoretical investigation was made of tandem hydrofoil 
arrangements, based on the lifting-line theory. The conclu- 
sions which follow apply to only the longitudinal behavior, 
inasmuch as the computations made were insufficient to 
justify definite conclusions regarding the lateral motions. 

1. The longitudinal hydrofoil spacing should be as large 
as is feasible. 

2. The rate of change in lift with change in depth of im- 
mersion of the hydrofoils should be large. Dihedral appears 
to be advantageous, if the hydrofoil is partly immersed, be- 
cause with dihedral there is a larger rate of change of lift 
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(a) r=30°; S 2 =Si; *=0; /=20.0ci; z„=1.74ci; Kr= 6.67ci; zi = 5.00ci; x\= 7.00ci. (b) l’=30°; S 2 =Si; e=0; I = 100.0ci; z. = l.74ci; NV=6.67ci; zi = 5.00ci; zi=35.00ci. 

Figure 17 Indictel responses for unit Cz and unit C m disturbances. 





Figure 18.— Significance of longitudinal hydrofoil spacing on a typical flying boat. 


with change in immersion. The rate of change of lift with 
immersion will be insufficient for hydrofoils with no dihedral 
unless the area is composed of several panels in a multiplane 
arrangement. 

3. The rear hydrofoil area should be as large as, or larger 


than, the front hydrofoil area if large variations in center-of- 
gravity location are to be accommodated when the longi- 
tudinal hydrofoil spacing is small (of the order of 10 chords). 
With appreciably larger spacings, the arrangement with the 
main surface forward appears to be sufficiently stable and 
should be more efficient than the other arrangements. 

4. The choice of horizontal center-of-gravity location 
should be based on considerations of the resultant charac- 
teristics of the longitudinal motions and the hydrofoil load- 
ing. The location should not be ahead of the hydrodynamic 
center of the front hydrofoil, in order to avoid unde- 
sirable loading. The location should be as far ahead of the 
rear boundary of the stable region as is feasible without 
incurring objectionable oscillations. The, best compromise 
from this latter standpoint appears to be a location near the 
center of the stable region. For two equal hydrofoils in 
tandem the best location appears to be back about 35 percent 
of the distance between the hydrofoils. 

5. If the effects of power are neglected, the vertical center- 
of-gravity location appears to be of little importance, low 
locations being somewhat advantageous. 

6. A reduction in the pitching radius of gyration will cause 
an appreciable increase in the range of horizontal center-of- 
gravity location that will be stable. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., May 9. 194-7. 



APPENDIX 

METHODS OF THEORETICAL TREATMENT 


SYMBOLS 


X-, Y-, Z~ axes 


x, Y,z 

L, M, N 
7 /- axis 


7/ 

<>, e , 'P 

a, 13 


V 

a 

U, V, w 


V, ?, r 

W 

m 

kx, ky, kz 


Pic 

1 

2 

s 


s n 

c n 

bn 

A n 

r. 


rectilinear reference axes fixed in hydrofoil 
system, with origin located at center of 
gravity (The X-axis is alined in the direc- 
tion of the initially undisturbed motion. 
The initial position of the F-axis is di- 
. rected horizontally to the right. The 
Z-axis is directed downward.) 
forces along X-, Y-, and Z-axes, , respec- 
tively 

moments about X-, Y-, and Z-axes, respec- 
tively 

axis, directed vertically downward with 
respect to the earth from origin located at 
center of gravity of hydrofoil system 
displacement along Z'-axis 
angular displacements of reference axes from 
initial positions, radians (see fig. 19) 
angles, in radians, giving instantaneous 
orientation of reference axes with respect 
to path of motion (see fig. 19); thus a is 
angle of attack and [3 angle of sideslip at 
center of gravity 

linear velocity of center of gravity 
angular velocity of hydrofoil system about 
center of gravity, radians per second 
components of V along X-, F-, and Z-axes, 
respectively 

components of about X-, Y-, and Z-axes, 
respectively 

weight of hydrofoil system 
mass of hydrofoil system 
radii of gyration of hydrofoil system about 
respective reference axes 
density of water 

subscript used to designate front hydrofoil 
in a system of two hydrofoils in tandem 
subscript used to designate rear hydrofoil in 
a system of two hydrofoils in tandem 
total projected area of immersed part of 
hydrofoil system under conditions of 
steady undisturbed motion 
total projected area of wth hydrofoil 
chord of nth hydrofoil 
span of nth hydrofoil 
aspect ratio of nth hydrofoil 
dihedral angle of nth hydrofoil, in radians 
unless specified otherwise 


r 




e 



dihedral angle when angle is same for all 
hydrofoils in system 

angle of attack at hydrodynamic center of 
nth hydrofoil, radians 
induced angle of attack at hydrodynamic 
center of front hydrofoil, radians 
downwash angle at hydrodynamic center of 
rear hydrofoil, radians 
rate of change of e with a 
rate of change of e with qc x jV 
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angle of sideslip at hydrodynamic center of 
nth hydrofoil, radians 
nondimensional rolling velocity at hydro- 
dynamic center of nth hydrofoil, based on 
local rolling velocity in radians per 
second, b n , and V 

nondimensional yawing velocity, with defi- 


nition similar to that for 



lift on hydrofoil system, measured at 
center of gravity in direction perpendicu- 
lar to V and converted to coefficient form 


by dividing by k p w V 2 S 


lift on nth hydrofoil, measured at hydro- 
dynamic center of hydrofoil under con- 
sideration in a direction parallel to C L 
and converted to coefficient form by 

dividing by^- p a V 2 S 


lift on nth hydrofoil, measured at hydro- 
dynamic center of hydrofoil under con- 
sideration in direction perpendicular to 
local relative motion and converted to 

coefficient form by dividing by p u V*S n 


drag on nth hydrofoil, measured at hydro- 
dynamic center of hydrofoil under con 
sideration in direction parallel to local 
relative motion and converted to coeffi- 
cient form by dividing by 2 p w V 2 S n 


weight of hydrofoil system converted to 
coefficient form by dividing by ^ p m V 2 S 


side force on hydrofoil system, measured a( 
center of gravity in direction of Y-axis 
and converted to coefficient form by 

dividing by ^ p,„V 2 S 


side force on nth hydrofoil, measured at 
hydrodynamic center of hydrofoil under 
consideration in direction parallel to 
Y-axis and converted to coefficient form 

by dividing by ~ p w V 2 S n 


coefficient of Z-force, with definition similar 
to that for C Y 

rolling moment about A'-axis, converted to 
coefficient form by dividing by ^ p w V 2 Sbj 

rolling moment at hydrodynamic center of 
nth hydrofoil about axis parallel to 
X-axis, converted to coefficient form by 

dividing by ~ p m V 2 S n b n 


C m 

G n 

(cy„ 

Xl 


■To 


l 


Zn 


Zo 


n 


Pc 


Pi, 

Ky 

Kx, Ky, 


t 


s 


c 


S b 

X 


pitching moment about Y-axis, converted 
to coefficient form by dividing by i p,„Sci 

coefficient of yawing moment, with defini- 
tion similar to that for Ci 
coefficient of yawing moment, with defini- 
tion similar to that for (C t ) n 

the derivative £><?£„/& ( y- 

A r -component of distance from center of 
gravity to hydrodynamic center of front 
hydrofoil, c i units 

A'-component of distance from hydro- . 
dynamic center of rear hydrofoil to center 
of gravity, Ci units 

distance between hydrodynamic centers of 
the two hydrofoils measured parallel to 
A r -axis, Ci units 

Z-componcnt of distance from center of 
gravity to hydrodynamic center of ?rth 
hydrofoil, c ; units 

operating depth; distance from water sur- 
face to hydrodynamic center of nth 
hydrofoil during steady undisturbed mo- 
tion, c, t units 

operating depth when depth is same for all 
hydrofoils in system 

parameter of /till hydrofoil used to deter- 
mine value of b(Cy) n /c)(^\ 

Y-componcnt of distance from hydro- 
dynamic center to centroid of lift on one 
panel of nth hydrofoil, b n units 
vertical displacement of center of gravity 
during disturbed motions, c t units 
vertical displacement of hydrodynamic 
center of nth hydrofoil during disturbed 
motions, c n units 

mass of hydrofoil system, 2 p a Sci units 

mass of hydrofoil system, ^ p w Sbi units 

radius of gyration about Y-axis, c i units 
radii of gyration about X- and Z-axes, 
respectively, hi units 
time, seconds 

time, c\jV units (To convert nondimen- 
sional time into second units use t 

The s c time scale may alternatively be 
converted into distance traversed if 
values of s c arc multilpied by Ci.) 
time, bi/V units 

stability root, with various numerical sub- 
scripts used to distinguish the different 
roots 
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m 


C z (s r ) 


k t 

k 2 , k-i 


disturbance function; a Z-force of variable 
magnitude, time history of which is indi- 
cated by form of function (The complete 
description of any arbitrary disturbance 
acting on the hydrofoil system may be 
expressed by use of this and the additional 
disturbance functions M(t), Y(t), L{t), 
and N(t), with definitions similar to that 
for Z(t).) 

nondimensional disturbance function, simi- 
lar to Z(t) but with force expressed in 
coefficient form and with time in nondi- 
mensional units (Similar definitions ap- 
ply to C m {s c ), Cy (s 6 ), Ci(s b ), and C n (s,,).) 

indicial responses giving motions a, z ' , and 
0, respectively, caused by sudden applica- 
tion of unit C z disturbance to hydrofoil 
system 

indicial responses giving motions in a, z' , 
and 6, respectively, caused by sudden 
application of unit C m disturbance to 
hydrofoil system 

empirical constant used to determine 
value of b(C L )„/bz' n 

empirical constants used to determine value 


of b(C L )„/ba n 

k it ha empirical constants used to determine value 

oib(C D )Jba n 

The abbreviations li.c. and c.g. are used herein for hydro- 
dynamic center and center of gravity of hydrofoil system, 
respectively. 


LONGITUDINAL EQUATIONS OF MOTION 


The longitudinal motions of the hydrofoil system arc 
referred to the system of axes described in the list of symbols. 
The choice of axes that correspond to those customarily 
employed in studies of airplane stability should facilitate 
extension of the present hydrofoil theory to include the 
effects of aerodynamic surfaces. The equations of motion 
are based on the assumption that the hydrofoil system can 
be replaced by a particle at its center of gravity having a 
mass m and radii of gyration k x , k Y , k z about the respective 
reference axes equal to those of the hydrofoil system. The 
analysis is also based on the assumption that the velocities 
V in the direction of motion and u along the X-axis are 
constant and that departures from the initial conditions of 
motion arc small. The further assumption is made that the 
longitudinal displacements 7', 6, and along the Z-axis, in 
the plane of symmetry of the hydrofoil -system,, are inde- 
pendent of the lateral motions involving the displacements 
<t>, ip, and along the Y-axis. This assumption yields satis- 
factory theoretical predictions of the motions of airplanes in 
normal flight and appears warranted, based on the nature of, 
the deviations involved, in the treatment of the longitudinal 
motion of hydrofoils. Its application to the lateral motions 
of hydrofoils is made with reservations, as mentioned in the 
main text. 


By the use of D’Alemberts principle, the following 
equations of equilibrium at the center of gravity are written 
for the forces and moments involved in the longitudinal 
motions: 


S (1) 


where Z(t) and M(t) are arbitrary disturbance functions. 
The equations have the same form as the familiar equations 
of longitudinal motion for an airplane, except for the 
addition of derivatives with respect to Z' and 9. The 
equation of equilibrium involving the X-force is omitted 
because u is assumed constant. Equations (1) can be 

. , . dZ w . dd , 

simplified by using w=-^> a ~y’ an d <Z = ^ to £ lve 


m 


„da r/ d9 bZ ry, bZ bZ ddbZ 

V dt~ mV dt =a ba +z b^ d be+dibi +Z{t) 

y ( 2 ) 


, „ d 2 8 bM bM bM.dd bM . , , , . \ 

m ky- d(2 =a da +Z dZ? +0 d0 + df J 


If equations (2) are rewritten in a nondimensional form, the 
solutions obtained will be general in character. The 
method used to make the various terms of the equations 
nondimensional involves expressing all angles in radians, 
all forces and moments in the standard NACA coefficient 
forms 


C z = 


c =- 

I 


z 

P,Y 2 S 

(3) 

M 

(4) 


; p u y 2 sc 1 


all lengths in terms of the chord Ci of the front hydrofoil, 
all times in terms of the time Ci/1 7 required for the system 
to traverse the distance c, along the path of motion, and 

the mass in terms of pu.Sc i units. The nondimensional 

quantities of mass p c , time s c , vertical displacement z', and 
radius of gyration K y about the Y-axis thus bear the follow- 
ing relations to the corresponding dimensional quantities: 


m 

Pc— 1 

*2 Pwk>C\ 

(5) 

t 

8e ~«i/V 

(6) 

s' = z - 

(7) 

Cl 

K •= — 
5 c. 

(S) 
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In equations (3) to (8), p w is the density of water and S is ' 
the total projected hydrofoil area in the hydrofoil system. 
The nondimensional form of equations (2) becomes 


M 


da 

dd\ 

<>Cz 

1 / dCz . a 

dC z 

ds c 

dsj 

= a 

Oa 

f2 d7 +9 

dd 

p c K\ 

0 d 2 6 

dC m . 


dC m 

r ds<? 

— a~< 1 

da 

■ z +e 

de 


dd ' ds c ^ gCi 


° V 




■CM 


Also, from geometric considerations, 

dz' 


ds, 


= a — 6 


( 10 ) 


In equations (9), C z (s c ) and C m (s c ) are functions of non- 
dimensional time that describe the application of disturbing 
force and moment coefficients to the hydrofoil system. The 
methods used to make the terms of equations (9) non- 
dimensional have the advantage that the nondimensional 
equations obtained retain the same form as the original 
force equations, consequently the physical significance of 
the nondimensional equations should be more readily evident. 
Solutions of motion obtained from equations (9) are likewise 
nondimensional and may be considered as proportions, 
applicable to all similar hydrofoil systems, and capable of 
conversion to customary engineering units in any given 
case by use of the characteristic dimensions c i and V perti- 
nent to the specific design. 

Stable regions and stability roots for the longitudinal 
motions can be obtained from equations (9) in conjunction 
with equation (10) by methods discussed in reference 4. 
The stability equation for the longitudinal motions has the 
form 

aD i + bD 3 + cD 2 +dD+e = 0 (11) 


Boundaries for the stable regions were obtained from the 
conditions 

( bc—ad)d — b 2 e — 0 (12) 


for the oscillatory boundary and 


e=0 


(13) 


for the divergence boundary. The quantities involved in 
equations (12) and (13) are the coefficients of equation (11), 
which in turn are functions of the factors of equations (9) 
and (10). Thus, 


fi _dCz fdCz,dC£\ 

dz' \ dd da J dz' \dd da J 


(14) 


Equation (12) is the familiar Routh’s discriminant, but its 
expression in terms of the factors in equations (9) and (10) 
is considered too lengthy to be presented here. 


LONGITUDINAL DERIVATIVES 

Values must be assigned to the various partial derivatives 
appearing in equations (9) before the equations can be 
solved. No experimental values for the derivatives were 
available; hence computed values were used. The com- 
puted derivatives were evaluated on the basis of experimental 
hydrofoil data obtained from results of tests made in the 
Langley tank no. 1 at various immersions and speeds. A 
discussion of the methods used to compute the various 
derivatives follows. Data presented in connection with 
the discussion are for hydrofoils of rectangular plan form 
and tips, with an aspect ratio of 6, and operating at a veloc- 
ity of 20 feet per second. Experimental results indicate 
that, for a given angle of attack, marked changes in the lift 
and drag coefficients of hydrofoils occur with changes in 
speed. The values of the derivatives would undoubtedly 
be equally affected by any pronounced change in speed 
from that assumed in the investigation. 

Change in Z-force with vertical displacement of the 
center of gravity dC z /dz '. — If the center of gravity moves 
downward, the hydrofoils are immersed deeper in the water. 
Experimental results indicate that an increase in the depth 
of immersion of a hydrofoil is accompanied by an increase in 
the magnitude of the lift obtained. The increase in lift is 
proportional to, and of the same sign as, the initial lift. 
Thus, 

(15) 

Values of k t are given in figure 20 (a) for a dihedral angle of 
0° and in figure 20 (b) for dihedral angles of 20° and 30°. 
The value of /e, depends on the normal operating depth z„ n 
of the hydrofoil. The discontinuities in the curves of 
figure 20 (b) coincide with the point where the tips of the 
hydrofoil break the surface. In figure 20 (a) and subse- 
quent figures, (C L ) n is based on the total area of the hydro- 
foil instead of the immersed area and z 0ji is measured in 
chord lengths of the particular hydrofoil under consideration. 

The value of dC z /dz' for a complete hydrofoil system is 
the negative sum of the values of dC L Jdz' for the individual 
hydrofoils. The values of dC L Jdz' for the various hydro- 
foils are derived from the d(C L )„/dz' n values obtained from 
figure 20 (a) or 20 (b) by making proper allowance for the 
different areas and chords that are used to make the various 
terms nondimensional. 

Change in Z-force with angle of attack dCz/da. — The value 
of the derivative dC z /da is the negative sum of the values of 
dC L Jda (that is, the slopes of the lift curves) for the individ- 
ual hydrofoils. As in the case of dC z /dz' , differences in the 
areas used in forming the coefficients must be taken into 
account when the addition is made. The slope of the lift 
curve depends on the depth of immersion of the hydrofoil. 
Typical variations of the slope are given in figure 21 (a) for 
0° dihedral angle and in figure 21 (b) for various dihedral 
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(a) -<4. n =6; NACA 66, S-209 section; r n =0°; c„=5 inches; V=20 feet per second. 

(b) /1„=6; NACA 16-509 section; r n =20° or 30°; c n — 5 inches; V=20 feet per second. 

Figure 20. — Variation of k\ with z on for hydrofoil having rectangular plan form and rectangular 

tips. 



(a) .An=6; NACA 66, S-209 section; r„— 0°; c n =5 inches; T r =20 feet per second. 

(b) >ln=6; NACA 16-509 section; r„=0°, 20°, or 30°; c„=5 inches; V=20 feet per second. 

Figure 21.— Variation of d(CL)nldx n with z on for a hydrofoil having rectangular plan form 
and rectangular tips. 


In general the value of a 2 is less than that of a (measured at 
the center of gravity) by the amount of the downwash angle e 
at the rear hydrofoil. The slope of the lift curve for the 
rear hydrofoil must be corrected for downwash to give the 
required slope with respect to a. The correction is applied 
by multiplying the slope obtained from figure 21 (a) or 21 (b) 
by the factor 1 — e a , where e a has some value in the range 




2 bjCi h 

Ax Qctx 


angles. When figure 21 is used to determine the slope of the 
lift curve for the rear hydrofoil, the value obtained is with 
respect to the local angle of attack <x 2 at the rear hydrofoil. 


( 16 ) 
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In equation (16), A x is the aspect ratio of the front hydrofoil 
and &(&)i/&ai is the lift-curve slope obtained from figure 
21 (a) or 21 (b) for the front hydrofoil. 

Change in Z-force with pitch attitude bC z /bd.- — A change 
in the pitch attitude of the hydrofoil system will cause a 
differential change in the depth of immersion of the hydro- 
foils. The effect on the Z-force may be estimated from the 
geometry of the system and the data of figure 20; thus, for 
two tandem hydrofoils 


bC z _S l d(C L ) s S 2 c, d(C 7 ,) 2 
be ~ S bz\ Xl S c 2 dz' 2 


(17) 


Change in Z-force with pitching velocity bC z /b ^ — The 


chief effect of a pitching velocity about the center of gravity 
of the hydrofoil system is to cause a change in local angle of 
attack at each hydrofoil. The change in effective camber for 
the pitching hydrofoil introduces a small additional compo- 
nent of vertical force. (See reference 5.) The total effect for 
two hydrofoils in tandem may be assumed to be 


where 






gi &(&), 

S da. 


(— Xi + 0.5) 


(18) 


(19) 


< \=§^(* 1 -'' +0 - 5 ir:) < 20 > 

In equation (20), b{C L ) 2 jba 2 is the lift-curve slope for the 
rear hydrofoil, based on the local angle of attack a 2 ; x l and 
x -2 are the AT-corn ponents of the locations of the front and 
rear hydrofoil hydrodynamic centers from the center of 
gravity expressed in terms of c x ; and e„ indicates the rate 
of change of downwash angle at the rear hydrofoil with 
change in nondimensional pitching velocity qcJV. The 
value of e, will be in the range 


2 C L 


( 21 ) 


Change in pitching moment with vertical displacement of 
the center of gravity bC m /bz'. — The changes in lift, mentioned 
in the discussion of the change in Z-force with vertical 
displacement of the center of gravity, produce moment 
changes about the center of gravity, the magnitude of which 
depend on the JV-components of the distances of the hydrofoil 
hydrodynamic centers from the center of gravity. The drag 
also increases with deeper immersion of the hydrofoils. 
Analysis of the data obtained in Langley tank no. 1 indicates 
that the change in drag can be expressed as 

*.(<?*) »*+*, ( 22 ) 

Values of k 2 and k 3 are given in figure 22 (a) for 0° dihedral 
angle and in figure 22 (b) for 30° dihedral angle. The drag 


changes multiplied by the Z-components of the distances 
from the center of gravity to the hydrofoil hydrodynamic 
centers give the drag contributions to the change in pitching 
moment. For two hydrofoils in tandem 

bC m _Sib(C L )i S 2 c 1 b(C L ) 2 S l b(C D ) l , S 2Cl b(C D ) 2 

bz'~S bz', Xl Sc 2 bz\ 2 S bz\ 2l S c 2 ds' 2 22 

(23) 


Change in pitching moment with angle of attack bC m /ba . — 
Physical considerations lead to the expression, for two 
hydrofoils. 


bC m _S l b(C L ) l S a b(C L ) 2 
b<x~ S bet! Xl S 1 ta) ba 2 



(i—o 



X 2 -\- 


d(Cz.) 2 

5 0(2 


[(ft), 

J 22 


b (Cp) i 

da, _ 


Zi 


(24) 


where ( Cd)i is the drag coefficient of the front hydrofoil based 
on the area of the front hydrofoil; (C D ) 2 is the drag co- 
efficient of the rear hydrofoil based on the area of the rear 
hydrofoil ; and z x and z 2 are the Z-components of the locations 
of the front- and rear-hydrofoil hydrodynamic centers from 
the center of gravity, expressed in terms of c,. 

The slope of the drag curve for each hydrofoil must be 
known to determine bC m /bct from equation (24). The 
empirical relation 

-^ )n =h(C r ,) n ~k, (25) 


was obtained from an analysis of the experimental data. 
Values of and k 3 'varied with the depth of immersion 
of the hydrofoils in the manner shown in figure 23 (a) for 0° 
dihedral angle and in figure 23 (b) for 30° dihedral angle. 

Change in pitching moment with pitch attitude bC m /bd . — 
The differential change in the depth of immersion of the 
hydrofoils introduced by a change in the pitch attitude of 
the hydrofoil system leads to variations in the lift and 
drag for each hydrofoil. These variations can be trans- 
lated into a variation in pitching moment about the center 
of gravity by use of the geometry of the hydrofoil system 
and equations (15) and (22). For two tandem hydrofoils 


bC m Si b(C L )i 2 S 2 ci b(C L ) 2 2 | 'S'i &(Cd) i _ 

be ~ S bz\ 1 S c 2 bz' 2 2 ^ S bz\ XlZ ‘ 


S 2 Ci b(C D )j 
S c 2 bz' 2 X2 " 2 


(26) 


Change in pitching moment with pitching velocity 

bC m fb The only important contribution to the pitching 

moment produced by a pitching velocity about the center 
of gravity is that associated with the change in lift on 
each hydrofoil as a result of the change in local angle of 
attack. Thus, 


b C m ___ S 1 b(C I ,) 1 
. qci S da, 

° V 


Xi 


S 2 . 

~g te — ««) 


biOh 

ba 2 X 2 


(27) 
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LATERAL EQUATIONS OP MOTION 

Equations expressing the equilibriums of the forces 
and moments involved in the lateral motions are written 
on the same assumptions as those used to obtain the longi- 
tudinal equations. The equations of lateral motion are 


d 2 Y 
m- 


i 2 Y . Tr ar. , roY „ 7 \, ar. ar lW ,' 

W+ mrV==V -^ + <p (j- + W)+p ^ +r^ + Y(t) 

7 2 d 2 <j> a L . ,bL . bL &L , r ... 

mk * r v bd +,f) bj+P bp +r -^ +L(t) 




bv 


bcj> F bp 


dr 


f ( 28 ) 


where Y(t), L(t), and N(t) are arbitrary disturbance func- 
tions. Equations (28) can be simplified by using v= 

„ y d<j> j dtp 
P= h> t )= d{’ and r= ~J7 to give 


dt 


-rrdp rrd'P O^Y 

mV Tt +mV Tt =fi bp + * 


(<>Y > w\+- d4> dF 4_ 1 

{bt^ w )+dtb£+ 


dip bY 
dt br 


+ Y(t) 


J 2 d 2 b bL bL dtp bL dip bL T 
mk - W^b^ +<t> bb + Tt bp + dtb7+ L{t) 

, 2 dV o 2>iV bN d<p bN dip bN 
mk/ W + *bb + diW + dt W +N ® 


(29) 


Equations (29) will next be written in a nondimensional 
form similar to that used for the longitudinal equations. 
Thus, all angles will be expressed in radians and all forces 
and moments in the standard NACA coefficient forms 


C--J 
- 2 P*V 2 S 

(30) 

w 

KJW — i 

2 PwV 2 S 

(31) 

c,= L 

2 p w V 2 Sb 1 

(32) 

r _ N 

n ‘ -j 

2 P,V 2 Sb x 

(33) 

Because of the different basis for forming the moment 
coefficients (cf. equation (4)) in the nondimensional lateral 
equations of motion, all lengths will be expressed in terms of 
the span of the front hydrofoil b x , all values of time in terms 
of the time b x /V required for the system to traverse the 


distance b x along the path of motion, and the mass in terms of 

~ p w Sbi units. The nondimensional mass p b , time s b , and 

radii of gyration K x and K z thus bear the following relations 
to the corresponding dimensional quantities: 




Py)Sb\ 


S b ~ 


bi/V 


K — kx 

Kx ~b, 

K — — 
Kz ~b x 


(34) 

(35) 

(36) 

(37) 


The nondimensional form of equations (29) becomes 


Pb 


\ds b ds b J 


w+ * 

dp bC Y 


I n \ I d<p bCy 
+ L w ) + dSb . p&r 
0 V 


'Cy(Si,) 


ds b . rbi 
° V 


K , <E$_o bCi , bCi d<(> bCi , dip bC, , 
x ds b 2 ^ d/3 b<p ' ds„ ‘pbi ds b . rb x 

0 V 

C t (s b ) 


jy 2 d^lp „ bGn # bCn , # bC n 

^ Kz ds b 2 ~P bp + </> b<p + ds b y^+38, rb + 

°v 0 V 


G n (s b ) 


y ( 38 ) 


where 0 Y (s b ), Ci(_s b ), and C n (s b ) are functions of nondimen- 
sional time that can be used to define the application of 
any lateral disturbance to the hydrofoil system. 


LATERAL DERIVATIVES 

In order to obtain a solution from equations (38), the 
various partial derivatives involved must be given numerical 
values. No experimentally determined values were availa- 
ble for any of the derivatives, and computed values were 
therefore used. Experience has shown that theoretical 
methods are unreliable for obtaining many of the lateral 
stability derivatives of airplanes. This fact, coupled with 
the additional complication of the presence of a free surface, 
suggests that theoretical computations of the derivatives 
for hydrofoils will be even less satisfactory. Elaborate 
theoretical analyses to obtain the values of the lateral 
stability derivatives of hydrofoils, therefore, appear to bo 
unjustified until experimental data are available for use in 
checking the accuracy of computed values. 
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For most of the lateral derivatives, the values of the deriva- 
tives were first computed with respect to the hydrodynamic 
center of the hydrofoil for motions at the hydrodynamic cen- 
ter; from the geometry of the hydrofoil system the deriva- 
tives at the center of gravity of the hydrofoil system for 
motions at the center of gravity were obtained. The follow- 
ing discussion will be mainly confined to methods of com- 
puting the lateral derivatives at the hydrodynamic center of 
the hydrofoil. Such derivatives can be readily converted to 
derivatives at the center of gravity of the hydrofoil system 
by the use of elementary mechanics when the geometry of 
the system is known. Numerical data presented in connec- 
tion with the discussion of the lateral derivatives were 
obtained from the same sources and the same operating condi- 
tions as those used in obtaining the longitudinal derivatives. 
The expressions derived are for the lateral derivatives of an 
“ideal” hydrofoil system without supporting struts. The 
presence of the supporting struts usually required will un- 
doubtedly have a large influence on the values of certain of 
the lateral derivatives. 

Change in F-force with sideslip dC y /d{3 . — During sideslip 
the effective angle of attack is differentially altered on each 
side of the hydrofoil, which changes the lift on each half in 
such a way that a component of side force is introduced. 
This effect is a function of the dihedral of the hydrofoil. In 
addition, the direction of the drag force is rotated to one side 
during sideslipping. The sum of these effects is 

3 r * tan r -- <c "> ■ <39) 


The lateral displacement of the centroid of lift from the 
juncture of the hydrofoil panels is given by the value of 
y Cji in figure 25. In order to- make y Cn nondimensional it is 
expressed in terms of twice the projected span of the banked 
panel. The new angle of attack of the panel after a change 
in bank is 

a=a 0 cos r 0 sec r (40) 


where the subscript 0 refers to the initial values for the hydro- 
foil panel, and r and a are the values of the dihedral angle 
and angle of attack of the panel after a change in bank. 
(Note that T = T 0 ±<t>, where the sign depends on whether 
the left or right panel is involved.) 

Equation (40) and the values of b(6i)„/c)a„ and y Cfl ob- 
tained from figure 25 can be used to determine the magni- 
tude and point of application of Cl for each banked panel. 
The value of C Y for the banked hydrofoil is then determined 
by rules of simple mechanics. The value of 3(7 y /3</> is 
obtained graphically by plotting the values of C y determined 
for several values of 4> and measuring the slope of the resulting 
curve. 

Change in F-force with rolling velocity bCyfd -y — An 


estimation of the value of the derivative 3CV/3 


pbi 

V 


was 


obtained on the assumption that the side force would be zero 
for rolling of the hydrofoil about its effective center of curva- 
ture in front elevation. The derivative for rolling about the 
center section of the hydrofoil can then be obtained by an 
expression of the form 


where (C Y )„ is the coefficient, based on S n , of the F-component 
of force at the hydrodynamic center of the nth hydrofoil 
and is the sideslip angle at the same point. The dihedral 
angle of the nth hydrofoil in radians is indicated by T n . 
The value of 3(CV,)n/3a„ required in equation (39) can be 
obtained from figure 21(a) or 21 (b), and the value of ( C D )n is 
given in figure 24(a) for 0° dihedral angle and in figure 24(b) 
for 30° dihedral angle. 

Change in F-force with angle of hank 3 CV/3<f>. — The value 
of the derivative 3Cy/3</> was estimated by treating each 
panel separately as a hydrofoil of which the dihedral angle, 
angle of attack, centroid of lift, lift-curve slope, and im- 
mersed area vary with angle of bank. The change in effec- 
tive aspect ratio, which should be small for small changes in 
bank angle, was neglected. The variation in dihedral angle 
and immersed area with angle of bank was obtained, by 
graphical methods, for banking about the center of gravity 
of the hydrofoil system. The changes in lift-curve slope and 
centroid of lift with dihedral angle were obtained from figure 
25. The value of d(C L ) n /da n in this figure is for a lift coeffi- 
cient based on the projected area of the hydrofoil while 
banked, rather than on the initial projected area, and with 
the lift measured vertically regardless of the bank attitude. 



3 (Cy) n 
3(3 n 


(41) 


The parameter r n is given in figure 26 for various dihedral 
angles. 

Change in F-force with yawing velocity dC Y /d y — The 
rb 

derivative 3t7 r /3 y was assumed to bo zero for yawing 

about the hydrodynamic center of the hydrofoil. 

Change in rolling moment with sideslip 36V 3/3. — The 
differential change in lift, produced on each panel of a hydro- 
foil during sideslip, introduces a component of rolling moment 
about the center section. An additional component of 
rolling moment arises because the point of application of the 
side force produced by sideslip lies above the center section. 
The sum of these effects is 


3 «?,)„ rf / 3 (<?,)„ 
fin ~ y ° n \ 3«~ 


r -+- l Sr tan r 


n ) 


(42) 


where y Cn is obtained from figure 25 and 3(C',.)„/3a, 1 from 
figure 21 (b). 
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(a) .<4n=G; NACA 66, S-209 section; r n =0°; c»=5 inches; V=20 feet per second. 

(b) /ln=6; NACA 16-609 section; I\»=30°; c„ = 5 inches; r=20 feet per second. 

Figure 24.— Variation of ( Cd)„ with ( Cc) n for a hydrofoil having rectangular plan form and 

rectangular tips. 


Change in rolling moment with angle of bank bCi/d<j >. — 
Increments of C L and C Y , caused by a change in angle of 
bank, can be computed by methods outlined in the discussion 
of dD y /d<£. These increments, when multiplied by appro- 
priate moment arms (expressed in span lengths), arc used to 
obtain a plot of Ci against 4>, from which the value of d(?i/d<£ 
is measured. 

Change in rolling moment with rolling velocity bCi d -y — 

Reference 6 gives —0.2 as an average value of the derivative 

dCz/d for a conventional airplane wing. The value for a 
I V 

hydrofoil will probably be somewhat smaller, but in the 
absence of experimental data the average value mentioned 
was used for rolling of the hydrofoil about its center section. 

/ 7'b 

Change in rolling moment with yawing velocity ~d(Ji r v' 
The average value 



was used for the derivative <iC, d - 


rb, 


V 


Reference 6 indicates 


that this value is suitable for wings with moderate taper, 
and the loss of lift on parts of a hydrofoil that approach the 
surface would result in a similar lift distribution if the hydro- 
foil had dihedral. 



0 0 10 20 30 40 50 

Hi, de 9 


Figure 25.— Variation of b{Ci.) n l<)a. n and ij cn with r„. 



Change in yawing moment with sideslip dC„/d(3 . — During 
sideslip the lift vector for each panel of a hydrofoil remains 
perpendicular to both the hydrofoil leading edge and the 
direction of motion. Hence, the projection of the lift 
vector on the horizontal plane rotates forward for the leading 
panel and rearward for the trailing panel. The resulting 
couple about the hydrodynamic center of the hydrofoil is 

-^-(Ca^tanr. (44) 

Change in yawing moment with angle of bank bC„/d<f >. — 
If, during banked motion of a hydrofoil, the centroid of drag 
for each panel is assumed to have the same location as the 
centroid of lift and if the additional assumption is made that 
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the variation of drag with lift is the same in the banked 
attitude as for zero bank, bC n /d(j> can be computed by 
methods similar to those used for c>CV/d<£ and d(7i/d<£. 

Change in yawing moment with rolling velocity ^ — 

The average value given in reference 6 for an elliptical dis- 
tribution of lift was used for the derivative bCJJd ?y- Thus 



(CW- 


d(Cx>) n 
(>CC n 


16 


( 45 ) 


The elliptical loading was assumed to approximate the loss 
in lift over the tip parts of a hydrofoil with dihedral and 
with the tips at the water surface. 

/ rb 

Change in yawing moment with yawing velocity d( 7 n /d y — 
The value 



(C!d)» 

8 


( 46 ) 


appears to be a suitable approximation to the expression 
given by Glauert for elliptical wings (see reference 6) and 
hence was used in the calculations. The selection of- ellip- 
tical loading was based on the same considerations as for the 

derivative c )C n jd^- 
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Positive directions of axes and angles (forces and moments) are shown by arrows 
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bol 
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0 

v 
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z 
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xL 
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Absolute coefficients of moment 

ri p ^ p N 

^ l ~qbS Vm ~qcS Ln -$bS 

(rolling) (pitching) (yawing) 


Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D Diameter 

p Geometric pitch 

p/D Pitch ratio 

V' Inflow velocity 

V, Slipstream velocity 

T Thrust, absolute coefficient C T — 

Q Torque, absolute coefficient C Q = 



P 


G, 

T 

V 

n 

pn 3 D* 


Q 


pn 1 !? 



Power, absolute coefficient C P = 


Speed-power coefficient = 




pn 3 Z ? 5 


pYI 

Pn 2 


Efficiency 

Revolutions per second, rps 
Effective helix angle ==tan _1 ^ 2 ~^ 


5. NUMERICAL RELATIONS 


1 hp=76.04 kg-m/s= 
1 metric horsepower= 
1 mph=0.4470 mps 
1 mps=2.2369 mph 


550 ft-lb/sec 
=0.9863 hp 


1 lb=0.4536 kg 
1 kg=2.2046 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 



>r-i 




